Quantum Hall effects in fast rotating Fermi gases with anisotropic dipolar interaction 
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We investigate fast rotating quasi-two-dimensional dipolar Fermi gases in the quantum Hall 
regime. By tuning the direction of the dipole moments with respect to the z-axis, the dipole- 
dipole interaction becomes anisotropic in the x-y plane. For a soft confining potential we find that, 
as we tilt the angle of the dipole moments, the system evolves from a, u — 1/3 Laughlin state with 
dipoles being polarized along the z axis to a series of ground states characterized by distinct mean 
total angular momentum, and finally to an anisotropic integer quantum Hall state. During the 
transition from the fractional regime to the integer regime, we find that the density profile of the 
system exhibits crystal-like structures. We map out the ground states as a function of the tilt angle 
and the confining potential, revealing the competition of the isotropic confining potential and both 
the isotropic and anisotropic components of the dipole-dipole interaction. 

PACS numbers: 03.75.Ss, 73.43.-f, 73.43.Nq 



I. INTRODUCTION 



The quest for quantum computer with intrinsic fault 
tolerance spurs recent interest in searching for exotic 
fractional quantum Hall (FQH) states that support non- 
Abelian anyons [1]. While it is easy to write down a trial 
wave function with highly nontrivial statistics, the real- 
ization of it in two-dimensional electron gases (2DEGs) is 
not simple. Apart from the technical difficulties of sam- 
ple preparation and low operating temperature, the lack 
of effective control on the interaction between particles is 
a major concern. In a realistic 2DEG, electrons interact 
via long-range Coulomb interaction, which may be modi- 
fied due to the presence of an adjacent gate in the case of 
graphene or the effect of Landau level (LL) mixing, which 
introduces effective three-body interaction (among oth- 
ers). In particular, recent theoretical and experimental 
studies [2-13 suggest that the perturbative modification 
of the interparticle interaction can have significant ef- 
fects on the stability of FQH states. Therefore, the ques- 
tion of how topological order evolves with interparticle 
interaction remains an interesting question with growing 
experimental capabilities of controlling microscopic pa- 
rameters. 

Ultracold atomic gases provide an ideal platform for 
simulating quantum many-body systems [6[. The re- 
alizations of FQH states in ultracold Fermi gases have 
been discussed in the presence of, for example, a rapidly 
rotating trap 0, |1] or a laser-induced geometric gauge 
field For identical fermions, s-wave interactions van- 
ish due to the Pauli exclusion principle. Unless in the res- 
onance regime, p-wave interactions in a single-component 
Fermi gas are typical very small. Nevertheless, signifi- 
cant interactions can still be introduced by using atoms 
or molecules with strong dipole-dipole interactions [lol - 
[12} . The FQH effects in a two-dimensional (2D) dipolar 



Fermi gas with isotropic dipole-dipole interaction have 
been studied in Refs. [l^, The system has been 

shown to undergo transitions from an integer quantum 
Hall (IQH) state to a = 1/3 Laughlin state, and to a 
Wigner- crystal state by increasing the rotation frequency. 
However, we are not aware of any work on the FQH ef- 
fects in the presence of anisotropic interaction, when the 
dipoles are not oriented along the rotation axis. 

In the present work, we study the FQH effects in a fast 
rotating quasi-2D gas of polarized fcrmionic dipoles. By 
tilting the direction of the dipole moments with respect 
to the rotation axis, we can tune the dipole-dipole inter- 
action to be anisotropic on the plane of motion. Starting 
from a Laughlin state with isotropic dipolar interaction, 
we investigate the ground state properties by varying 
the tilt angle of the dipole moments. We find that for 
small tilt angle the ground state can be approximately 
described by a FQH state. However, as one further in- 
creases the tilt angle, the ground state deviates from the 
FQH state significantly such that a crystal-like pattern 
emerges in the density profile of the gas. The ground 
state of the system eventually becomes an IQH state 
when dipole moments are aligned in the 2D plane. For 
a soft confining potential, the IQH state is noticeably 
anisotropic. We map out the phase diagram in the pa- 
rameter space spanned by the tilt angle and the strength 
of the confining potential. The results can be explained 
by the competition of the isotropic confining potential 
and both the isotropic and anisotropic components of 
the dipole-dipole interaction. 

This paper is organized as follows. In Sec. [Hi we 
present our model and calculation for relevant matrix 
elements of the model Hamiltonian. Section IIIII briefly 
covers the FQH states with isotropic dipolar interac- 
tion for later comparisons. In Sec IIVI we investigate 
the ground state structure in the presence of anisotropic 
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dipole-dipole interaction in a weak confining potential. 
The full phase diagram is presented in Sec. [V] We con- 
clude our discussions in Sec. IVll 



of the axial harmonic oscillator, with a wave function 



= 7r-l/4„-l/2„-2:V(2g") 



, where q ~ \J h/{^u)z)- In- 



tegrating out the variable z from Eq. ([T]), we obtain the 
Hamiltonian for the quasi-2D system as 



II. MODEL 

We consider a system of N spin polarized fermionic 
dipoles trapped in an axially symmetric potential 

f/(r) = i/i(L.V+t^V + ^'^'), 

where ^ is the mass of the particle, uj and ujz are the 
radial and axial trap frequencies, respectively. The trap- 
ping potential rotates rapidly around the z-axis with an 
angular frequency VI {< i^)- We further assume that the 
dipole moments d of all particles are polarized by an ex- 
ternal orienting field which is at an angle 9 about the 
z-axis. Since the s-wave coUisional interaction vanishes 
for spin polarized fermions, particles only interact with 
each other via a dipole-dipole interaction. If the orient- 
ing field corotates with the trapping potential, the dipo- 
lar interaction becomes time-independent in the rotating 
frame, i.e., 

V(r)=c,l/i^°)(r), 

where Cd — / (47reo) or j (47r) for, respectively, elec- 
tric or magnetic dipoles, with £o (mo) being the vacuum 
permittivity (permeability). The spatial dependence of 
V(r) can be described by 



v. 
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(r) = — P - 3(z cosd + x sin df\ 



Here, without loss of generality, we have assumed that 
the dipole moments are polarized in the x-z plane of the 
rotating frame. We can tune the dipolar interaction by 

introducing a tilt angle Q such that vf"^^ (r) is isotropic 
(anisotropic) on x-yplane for 6* = (6* ^ 0). In the 
rotating frame, the Hamiltonian of the system becomes 



3D 



E 



2^1 



u{vi) - nil 



i<j 



(1) 



where = xpy — ypx is the z component of the orbital 
angular momentum. 

Under the condition ^ w, the system can be re- 
garded as quasi-2D. As a result, the motion of all par- 
ticles along the z-axis is frozen to the ground state 



H2B - 



2^ 



(2) 



where p = (a;,y), is the unit vector along the z-axis 
and 
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(27rg2)i/2 



dze-^^/(2^^)l^i^^)(p,z). (3) 



The first term on the righthand side of Eq. ^ repre- 
sents the single-particle Fock-Darwin Hamiltonian in the 
symmetric gauge p^[l6| , which can be solved exactly to 
yield eigenenergies [8| 



h{Lj — r2)n_|_ + h{uj + r2)n_ -(- fvjj, 



(4) 



known as the Fock-Darwin levels, where the quantum 
numbers and n_ are two non-negative integers. In the 
fast rotating limit i7 — > w, the Fock-Darwin levels mimic 
the LLs with a level spacing 2MI. Throughout this work, 
we assume that the interaction energy is much smaller 
than the LL spacing, such that particles only occupy the 
highly degenerate lowest Landau level (LLL). 

To proceed further, it is convenient to introduce a set 
of dinicnsionless units: h for angular momentum, t — 
yjhj {2[Ujj) for length, and Cdjl^ for energy. The wave 
function of the LLL can then be expressed as 



V'm(p, V') 



pm^imip^-p /4 

V27r2"m! 



(m > 0), 



describing a state with an angular momentum mh. 
Within the LLL formalism, the Hamiltonian ([2]) in the 
second quantization reads 



H = aL^ 



Vi234(6')/li/rL2/rn4/i 



(5) 



mim2Tn3m4 



where /^^ is the fermion creation operator that cre- 
ates a particle in state ipm and = ''^fmfm the 
total angular momentum. The dimensionless quantity 
a — h{uj — VL)i^ /cd characterizes the relative strength of 
confining potential with respect to interaction. In the 
presence of anisotropy, the interaction matrix elements 
are 
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3 cos2 e - 1 



234 



J1234 ~ 4/C1234 Smi+m2,m3+m4, + Sin 6'/Cl234<5mi +m2 ,m3+m4±2 



(6) 



where 



"^1234 — 



■|m3 — mi| — |m4 — m2 



7<.l'frr7<ll 



2%/27r(j 2(l™3-™il+l™'i-'"2|)/2 W [m{'3]![TO5"4]!' 



dipole moment. In the following sections, we will explore 
the quantum states of the system in the parameter space 
{N,9,a). Our main focus is on the parameter ranges 
of TV < 10, < 6* < |, and 0.01 < a < 0.1. Unless 
otherwise stated, the value of q is chosen to be 0.01£. 



>7l234 = / dpp 



I m 3 — ni 1 1 + 1 '"4 - »Ti 2 I + 1 „ - P 



1234 



, . 7- I"i3-mi I / P \ T\mi-in2\( P i 

'yj ™?4 ly'' 



/ ■ lm3-mi| + |m4-m2|+2 -p^+g^pV2 



^*-"f^)^:;.r'"'(T)"K^ 



where mf^ — minimi, rrij) and = max(mi, mj). 

L^(-) is the associated Laguerre polynomial and erfc(-) 
is the complementary error function. Clearly, for 9 — 
Vi234 is nonzero only when mi + m2 — — 1714 — 0, 
indicating that the total angular momentum Lz is con- 
served in the isotropic case. However, when the interac- 
tion becomes anisotropic (0^0), V1234 are nonzero when 
mi -|- m2 — — = or ±2. 

Hamiltonian ([S]) contains three parameters: the total 
number of particles N = J2m fmfm, the relative strength 
of the confining potential a, and the tilt angle 9 of the 



III. QUANTUM HALL STATES WITH 
ISOTROPIC DIPOLAR INTERACTION 

Let us first assume that the dipolar interaction is 
isotropic in the x-y plane, which corresponds to 9 = 
in the Hamiltonian ([5]). In this case, the total angular 
momentum is conserved. Therefore, one may numeri- 
cally diagonalize the Hamiltonian ^ in the subspace of 
a given total angular momentum M to obtain 



H 



(N) 
M.n 



= E 



(N) 
M,n 



(N) 
M,n 



(7) 



where E 



(AT) 



M,; 



and 



(JV) 



"Ain/ eigenenergies and eigen- 
states, respectively. The index n labels the state in the 
subspace of total angular momentum M with increas- 
ing eigenenergy, i.e., n — for the lowest energy state, 
n = 1 for the first excited states, etc. We emphasize that 
we present a numerically exact treatment of the dipolar 
interaction potential for a quasi-2D system with finite 
wave-function spread along the perpendicular direction, 
while the ideal 2D case has been studied previously [11] . 



TABLE I: Magic numbers in an isotropic TV = 10 system with various gs. The magic numbers obtained from the composite 
fermion theory are included for comparison. 
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In Fig. [U we plot the eigenenergy versus the total an- 
gular momentum for = 6 and a = 0. As a guide to 
the eyes, we have connected the lowest energy state in 
each total angular momentum subspace by a piecewise 
straight line, on which a series of shoulders are visible. 
The first state on the each shoulder, where a downward 
cusp appears in the spectrum, represents a possible can- 



didate for the global ground state of the system with 
increasing a. For a given a 7^ 0, only one of these states 
is the global ground state of the system; the correspond- 
ing total angular momentum is so called a magic number. 
The properties of these states have been studied exten- 
sively. Laughlin first noticed that the lowest-energy state 
of the M = 3TV(iV - l)/2 is closely related to the frac- 
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FIG. 1: (Color online) Energy spectrum of Hamiltonian ((Sj 
with A'' = 6, Of = 0, and = 0. The solid (red) line connects 
the lowest energy state in each total angular momentum sub- 
space as a guide to the eyes. The smallest angular momentum 
state on the each shoulder of the line may become the global 
ground state of the system as a varies. 
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FIG. 2: (Color online) Mean total angular momentum M 
versus the tilt angle 9 for a — 0.01 in systems with various 
sizes up to A'' = 8. For sufficiently large number of particles 
(A'' > 5), plateaus corresponding to distinct ground states 
develop on the curves. 



tional quantum Hall effect [17]. To translate these num- 
bers into asymptotic filling factors in the thermodynamic 
limit, Girvin and Jach proposed an explicit expression for 
these states in the finite system [18], 



Lo 
M' 



(8) 



where Lq = N{N — l)/2. In the context of quan- 
tum dots, Jain and Kawamura proposed an explana- 
tion of the magic numbers using the theory of compos- 
ite fermions 19, '20I, which has been further discussed 
in later works 21-23]. We list the magic numbers of 
our quasi-2D model for = 10 in Table U for several 
choices of q. For small q up to 0.1, our results are con- 
sistent with the results in 2D rotating Fermi gases with 
isotropic dipolar interaction studied earlier by Osterloh 
et al |13j . However, for larger q we observe a few discrep- 
ancies. Interestingly, the magic numbers 69 and 97, not 
showing for small q, are consistent with the prediction of 
the composite fermion theory. 



the angular momentum m, such that the diagonaliza- 
tion procedures are carried out in the space constructed 
from V'mS with ni < rricut- We have chosen a sufficiently 
large mcut such that our results presented in this work 
are not the choice of mcut- In numerical diagonalization, 
we take advantages of the fact that dipole-dipole inter- 
action only couples angular momentum subspaces that 
differ by an even angular momentum and diagonalize in 
the Hilbert space spanned by odd and even angular mo- 
mentum states separately. 

For a given set of parameters {N,a,d), the ground 
state wave function of the system is denoted by 
|4'(^)(a, ^^)) with the ground state energy £;(^)( a, u). 
Throughout this section, the strength of the confining 
potential is fixed at a = 0.01, such that the ground state 
is a I' = 1/3 Laughlin state for 6 — 0. The dependence 
of the results on a will be discussed in the next section. 
Due to the large size of the Hilbert space involved in 
the numerical diagonalization, We study systems of up 
to = 8 fermions. 



IV. QUANTUM HALL STATES WITH 
ANISOTROPIC INTERACTION 



Now we turn to the study of the quantum states in 
systems with anisotropic dipolar interaction. Since the 
total angular momentum is no longer conserved, one 
has to numerically diagonalize the Hamiltonian ([5]) in 
much larger Hilbert spaces. In practice, one may trun- 
cate the Hilbert space by introducing a cutoff, rricut, to 



A. State transitions induced by varying 9 

Even though the total angular momentum is no 
longer a good quantum number, its mean value. 



can be defined and we will see that it is sufficient to char- 
acterize the ground state of a system. In Fig. ^ we plot 
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the 6 dependence of M for systems with N < 8. We 
find that M always decreases with 0. As will be shown, 
this monotonically decreasing behavior of M{0) is be- 
cause the dipolar interaction becomes less repulsive as 9 
is increased, which reduces the size of the system [23 |. 
Interestingly, for sufficiently large number of particles 
{N > 5), plateaus develop on the M{9) curves. Roughly 
speaking, for 6 < 40° the mean total angular momen- 
tum decreases abruptly from one plateau to another as 
the tilt angle is increased, signaling sharp transitions be- 
tween ground states with distinct properties at various 9. 
For a concrete example, we examine in detail the transi- 
tions for the system of = 6. The sudden drops of M 
are observed as 



M : 45 



33.3° 



3Q Q° 

39.0 ^ 34.89, 



(9) 



where the numbers above the two arrows denote the criti- 
cal tilt angles. From the analysis presented in Sec. lIIH the 
first and second plateaus clearly correspond to the frac- 
tional quantum Hall states with filling factors = 1/3 
and 5/13, respectively, in the notation define in Eq. ^ 
by Girvin and Jach [l8| . The third plateau has an mean 
total angular momentum of 34.89, which is 0.3% smaller 
than 35, indicating that it mainly contains the FQH state 
with filling factor i/ = 3/7. As one further increases 9, M 
decreases smoothly toward an asymptotic value of 17 at 
9 = 90°. As shown in Fig. [51 similar features also appear 
in the systems with N = 5, 7 and 8. For A^ = 3 and 4, 
however, M always varies smoothly, consistent with the 
expectation that quantum phase transitions happen only 
in thermodynamically large systems; the absence of the 
sharp steps for A^ < 5 is the normal finite-size artifact. 

To gain insight into those transitions, we examine the 
energy spectrum of the system. In Fig. [31 we plot the 
low- lying energy levels versus the mean total angular mo- 
mentum for N = 6 for the isotropic and the anisotropic 
cases. In the isotropic case (9 = 0), the mean total angu- 
lar momentum M is simply the total angular momentum 
, which is a good quantum number. In both cases, 
energy levels group into clusters which are plotted with 
different colors. For convenience, we shall refer to each 
energy cluster using the angular momentum of the lowest 
energy level in the cluster. Even though the boundary 
between two adjacent clusters are not well defined for 
high-energy states, the lowest ones are clearly well sepa- 
rated. In Fig. ^a), we choose the range in which three 
energy-level clusters are visible. Among them, the low- 
est energy state in cluster-45 represents the ground state 
of the system for 9 = 0. This is the 6-particle Laugh- 
lin state as discussed in Sec. IIIII Apart from its mean 
total angular momentum, we can also observe the fea- 
ture of a Laughlin state also in the low-energy spectrum. 
First of all, the low-lying excitations are chiral, appear- 
ing only on the side of M > 45. In particular, at roughly 
M = 46 and 47, we find one and two states, respec- 
tively. They can be interpreted as the edge excitations of 
the ground state droplet, with their wave functions ap- 
proximated by the ground state wave function multiplied 
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FIG. 3: (Color online) Low-lying energy spectrum for an A'^ = 
6 system with a = 0.01 and a tilt angle of (a) 9 = 0° and (b) 
33.3°. Distinct colors are used to specify different energy- 
level clusters, which can be regarded as chiral excitations of 
the corresponding lowest-energy state. 



by a symmetric polynomial of the corresponding degree. 
Near M = 48, we expect three low- lying states but only 
find two; however, there is another level well above (near 
0.7), presumably due to the infiuence of the cutoff in the 
momentum space. The series of numbers of the low-lying 
states is consistent with the chiral Luttinger liquid theory 
and signifies the topological order of the corresponding 
ground state 



25| 



As one increases the tilt angle 9, all energy-level clus- 
ters move downward to lower energies because the dipo- 
lar interaction becomes less repulsive. Nevertheless, the 
counting of the low-lying excitations remains robust, sug- 
gesting the topological order is not destroyed by small 
anisotropy, as shown in Fig. [3jb). In addition, the clus- 
ters with lower angular momentum move faster than 
those with higher angular momentum, hence, for exam- 
ple, the lowest energy state in cluster-39 becomes the 
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FIG. 4: (Color online) Fidelity iPsg, as defined in Eq. ([10} , as 
a function of 6 for = 6 and a = 0.01. The inset magnifies 
the fidelity in the vicinity of the small dip around 9 — 49.3°. 



ground state of the system at = 33.3°. As one further 
increases 9, the lowest energy state in cluster-35 will be- 
come the ground state of the system. They have different 
low-energy excitation structure from that of the AI — 45 
ground state. In the case of the cluster- 39, in particu- 
lar, there are two chiral excitations at AM « 1, clearly 
different from the Laughlin case. 

The phase transitions induce by tuning anisotropy can 
be further confirmed by calculating the fidelity of the 
ground state wave function i27|] , 



(10) 



where 60 is a small quantity. The fidelity measures the 
similarity between two adjacent states in the parameter 
space 1271 . In the bulk of a single quantum phase, two 
states close in the parameter space have wave functions 
that are only perturbatively different, hence Tse is close 
to unity for sufficiently small 69. Near the phase bound- 
ary, two states that are close in the parameter space 
can have very different structure in their wave functions, 
therefore the fidelity can drop sharply in the quantum 
critical region, signaling a quantum phase transition in 
the thermodynamic limit. In Fig. |4l we plot the 9 depen- 
dence of J^se for a systems with N = 6 particles. Two 
transition points can be clearly identified and the critical 
9 values are consistent with those obtained in Fig. [21 A 
closer look at the Fse{9) (inset of Fig. [Jj) further reveals 
that there exists a third dip at 9 ~ 49.3°. The com- 
parison of the fidelity Tse for two different ^^s indicates 
that the similarity of the ground states decreases with 
the increasing distance between the parameter 59 near 
the dip. 
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FIG. 5: (Color onhne) Overlap integral Ofl^ [Eq. ([TT])] versus 
e for M = 45 (i^ = 1/3 Laughlin state), 39', 35, and 15 (IQH 
state) in the = 6 system with a = 0.01. 



B. Structure of the ground state wave function 
with anisotropic interaction 

To reveal the structure of the ground state wave func- 
tion [^''■^^(a, 6*)), let us calculate the overlap integral 



(11) 



between the isotropic and anisotropic ground states. 
Again, we present the data of the system with iV = 6 
particles. Figure [5] shows the dependence of on 
anisotropy for M = 45, 39, 35, and 15. As can be seen, 
the 0-axis is roughly divided into four regions, the bound- 
aries of which coincide with the three dips of the fidelity 
in Fig[4l 

In the first region 9 < 33.2°, the overlap of^^ is 
greater than 0.974, which indicates that the dominate 
contribution to the ground state wave function comes 



from the Laughlin state 



(6) 
45,0 



Nevertheless, close 



to the right boundary of this region, several states in 
A/ — 47 and 43 manifolds are mixed into the ground 
state wave function, such that of^ drops notably. The 
ground state wave function in region 33.3° < 9 < 39.6° 
mainly contains the states from M = 43,41,39, and 37 



manifolds. Particularly, the state 



(6) 
39,0 



provides the 



largest contribution to the ground state with 0.858 < 
Ofoio < 0.96. 

In the region 39.8° < 9 < 49.3°, the situation is much 
more complicated compared to those in the first two re- 
gions. At the left boundary, O^c^q is as high as 0.814, but 
it quickly drops to close to zero as the right boundary is 
approached. In fact, many states from M = 21 to 37 in 
the isotropic case contribute collectively to the ground 
state wave function. As a result, complicated structures 
are developed in the density profile of the system. In 
Fig. [SI we present four typical patterns in the density 
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(a) e = 40° 




(b) e = 43.2° 

O 



(c) e = 45.8° 

O 



(d) ^ = 49.3° 




FIG. 6: (Color online) Density profiles q(x, y) for various 6*8 
in tfie A'" = 6 system witli ct — 0.01. Tlie brightest points 
indicate the local majcima in the density of the system. 



profiles 

g{x,y) = ^(*W(a,^^)|/,U„,|vI/W(a,^^)^ 

mm' 

'XfPmix,y)i'rn'{x,y) (12) 

of tlie system in this region. For 9 = 40° [Fig. [B^a)], 
besides the commonly seen ring-shaped density, a vertical 
ridge appears along the y axis. Five local density maxima 
can be identified in the figure: 4 of them on the ring and 
the other at the center of the trap. Figure IH^b) shows 
the result for 9 = 43.2°, on which 6 local density maxima 
appear on the ring structure. Compared to the case at 
9 = 40°, the density profile is clearly stretched along 
the X axis, which also represents a generic trend for the 
density profile as 9 is increased. The reason behind this 
is because the interaction energy is lowered by stretching 
the gas along the direction of dipole moment. As 9 is 
further increased to 45.8° [Fig.lGjc)], the structure with 6 
local density maxima becomes more prominent, such that 
each of them almost becomes an isolated island, as in a 
crystal structure. In Fig.lH^d), the 6 density islands start 
to merge. We remark that these crystal-like structures 
appear as a result of the interference between the states 

'^M^ri^' system is struggling to maintain a balance 
between a large set of competing states. We note that 
our calculations are based on a finite size system, this 
region of competition may shrink in the thermodynamic 
limit, since more plateaus are developing as system size 
increases as shown in Fig. [21 

In the last region 49.3° < 9 < 90°, the most important 
contributions to the ground state wave function are pro- 



TABLE If: Largest 5 overlaps of tire anisotropic state 

(a, 90°)^ with the corresponding eigenstates in the 

isotropic case {9 — 0). The pair of (Af, n) labels the nth low- 
est eigenstate in the total angular momentum A/ subspace. 
The system has n — 6 particles and a = 0.01. 



{M,n) 


(15,0) 


(17,0) 


(19,2) 


(21,7) 


(23, 14) 




0.658 


0.608 


0.387 


0.196 


0.084 



vided by the states 



(6) 
15.0 



(6) 
17,0 



(6) 
19,2 



and 



(6) 
21,7 



IQH state 
that, at 9 



The probability of finding the system in either of the 
above 4 states is larger than 0.92 throughout this region. 
To understand the properties of the ground state in this 
region, we first consider the specific state at 6' = 90° . Ta- 

ble|TT]hsts the largest values of ojj*^^„(G.01, 90°). Clearly, 

the main components of [^'^^•'(0.01, 90°)) are the v = 1 

jg ^y and its edge states, which suggests 

90°, the ground state is an IQH state. To 
confirm this, we plot the density distribution q{x, y) of 
this state in Fig. [Tja). As can be seen, the surface den- 
sity for the elliptical plateau is exactly l/27r£^, which is 
identical to that of a = 1 IQH state. Since the fill- 
ing factor can be defined as v = 2Tr£^nf with n/ being 
the fermionic surface density [l^, the filling factor for 
the 9 = 90° state is thus 1. The homogeneous elliptical 
plateau of the density profile implies that our conclusion 
can be generalized to the thermodynamic limit. We plot 
the energy spectrum of the system in Fig. [Tj^b). It is 
well-known that the topological properties of the v ^ 1 
IQH state is labeled by the number of edge states, i.e., 
1, 1, 2, 3, 5, . . . for AL^ = 0, 1, 2, 3,4,.. ., which is exactly 
the case shown in Fig.[7jb), although in this case we have 
to replace AL^ by AAI. All the above evidences mount 
to the fact that the 9 = 90° state is a = 1 integer 
quantum Hall state. Since there is no phase transition 
observed in this region from various criteria, we conclude 
that the ground state for 9 > 49.3° can be characterized 
as an anisotropic integer quantum Hall state. 



C. Understanding the anisotropic dipolar 
interaction 

To understand how anisotropy in the dipolar interac- 
tion leads to the emergence of the anisotropic IQH state, 
we decompose the dipolar interaction into isotropic and 
anisotropic (on x-y plane) components as 



(3D) 



= K 



(3D) 



(3D) 



'7iso(6')- 



3z' 



(13) 



where r]iso{d) = (3cos2 6' - l)/2 and 77ani(6') = 3sin2 6'/2 
represent the strengths of the isotropic and anisotropic 
components, respectively. One should note that, to ob- 
tain Eq. (jl3|) . we have neglected the linear terms in z 
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FIG. 7: (Color online) Characterization of the anisotropic 
IQH state at e = 90° for = 6 and a = 0.01. (a) The den- 
sity profiles along the x- and y-axes highlight the anisotropic 
nature of the state. Near the center, the LLL is completely 
filled. The inset shows the ellipsoidal density profile on the 
x-y plane, (b) The low-lying energy levels can be counted as 
1, f , 2, 3, 5, . . . for AM = 0, 1, 2, 3, 4, . . ., respectively, which 
is consistent with the edge theory of an IQH state. 



as they vanish after we integrate out the variable z to 
obtain the 2D interaction potential. The properties of 
the system can be seen as determined by the compe- 
tition between vjj^^ and V^ ^ni'. Since ryiso (??ani) is a 
decreasing (increasing) function of the tilt angle, varying 
9 will change the relative strength of the isotropic and 
anisotropic components, which can give rise to different 
quantum phases. 

After introducing the decomposition Eq. (I13p . we ex- 
plore the contributions of the isotropic and anisotropic 
components of the dipolar interaction separately. To this 
end, we consider two fictitious systems, FS-1 and FS-11, 
in which the full dipolar interaction is replaced, respec- 
tively, by V'g jj^' and Vg^^^ . The corresponding Hamilto- 
nians take the same form as Eq. ^ except for that the 
interaction matrix elements are replaced by those calcu- 



lated using Vg^l^J or Vg^^^) . 

In FS-I, the strength of the isotropic interaction, riiso{9) 
decreases with 9. Therefore, increasing the tilt angle is ef- 
fectively equivalent to increasing the strength of the con- 
fining potential a for the full Hamiltonian with 9 — 0. We 
plot, in Fig. [SI the 9 dependence of the ground state an- 
gular momentum for system-I with N = 6 and a = 0.01. 
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FIG. 8; (Color online) Mean total angular momentum as a 
function of 6 for N = 6 and a = 0.01 in systems of various 
interactions. The solid, dashed, and dash-dotted lines corre- 
spond to, respectively, the real system with dipolar interac- 
tion, FS-I with the isotropic component only, and FS-II with 
the anisotropic component only. The vertical dotted line indi- 
cates the position of angle 9c, where the isotropic component 
becomes attractive. 



As 9 is varied, FS-I experiences the same transitions as 
those studied in Sec. lIIII For the first two transitions, the 
critical values of 9 roughly agree with those obtained us- 
ing full dipolar interaction potential. In particular, Vg ^^^^ 
vanishes at angle 9c = 54.74° and becomes attractive in 
the x-y plane for 9 > 9c- Consequently, the ground state 
becomes a. v = 1 IQH state for 9 > 9c- The reason 
that the transition to the IQH state occurs at a tilt angle 
smaller than 9c is due to the finite a used in Fig. [51 



In FS-II, we note that the system must be in the v = \ 
integer quantum state at = where Vg'^^J vanishes. As 
one increases 9, the mean total angular momentum M 
increases smoothly from 15 to roughly 21 (Fig. [S|), which 
suggests that the system always stays in the IQH state 
independent of the tilt angle, although the quantum Hall 
droplet is stretched gradually along the x axis. Further 
study on the ground state of FS-II can be carried out as 
those have been done in the previous subsection. 

From the above analysis, it becomes clear that the se- 
ries of ground state transitions induced by varying the 
tilt angle are mainly caused by the isotropic component 
of the dipolar interaction. The anisotropic component, 
on the other hand, changes the fine details of the ground 
state, as it mixes in states with different total angular 
momentum to the otherwise isotropic ground state, as 
is clearly exemplified in the anisotropic integer quantum 
Hall state. 
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^(degree) 

FIG. 9: (Color online) Mean total angular momentum 
M{0, a) of the ground state for a system with A'' = 6 par- 
ticles. The lower left corner is the z/ = 1/3 Laughlin phase, 
while the large (blue) region containing the upper right corner 
is the IQH phase, which crosses over to an anisotropic IQH 
phase at the lower right corner. 



V. GLOBAL PHASE DIAGRAM 

In Sees, mil and IIVI we have shown the ground 
state transitions induced by varying either the confining 
strength a or the tilt angle 9. We present the complete 
phase diagram in Fig. [3J in which we plot the mean to- 
tal angular momentum M as a function of 6 and a for 

= 6. The phase diagram is separated into regions with 
different value of AI, which are well defined along the a- 
axis, where the interaction is isotropic, as discussed in 
Sec. IIVI Within our numerical capabilities, we find the 
basic structure of Fig. [9] remains unchanged as system 
size N varies. 

We find that the Laughlin state with ly — 1/3 is robust 
for weak confinement and not too large tilt angle (< 30°), 
which assures that unintentionally introduced anisotropy 
in interparticle interaction is not important. On the other 
hand, FQH cannot survive in the large anisotropy of the 
dipolar interaction, when the isotropic component of the 
interaction becomes soft, as analyzed in Sec. IIVCI 

When the confining strength becomes stronger, the 
mean total angular momentum M becomes smaller, indi- 
cating the system of particles becomes denser and denser. 
The evolution of M is not smooth, but goes through a 
series of magic numbers, which can be interpreted by 
the corresponding filling factors. In the large confine- 
ment limit, the system develops into a maximum den- 
sity droplet with — 1. The state crosses over to an 
anisotropic maximum density droplet for large tile an- 
gles (at small confinement since the isotropic component 
of the interaction changes from repulsive to attractive in- 
teraction), as revealed in Fig. [71 which reflects the com- 



petition of the isotropic confining potential and both the 
isotropic and anisotropic components of the interaction. 



VI. CONCLUSION 

To summarize, we investigated the quantum Hall ef- 
fects in the LLL of a fast rotating quasi-2D Fermi gas 
with anisotropic dipolar interaction through exact nu- 
merical diagonalization. With the tilt angle of the dipole 
moment 9, we introduced a new control knob to explore 
the FQH effect in cold atomic systems. We studied in 
details the phase diagram of a finite-size system and con- 
cluded that phase transition is expected as the tilt angle 
6 varies in the thermodynamic limit. When the tilt angle 
is small, the ground state of the system can be described 
by a FQH state, whose filling fraction depends on the 
strength of the confinement and hence the average den- 
sity. At large tilt angle, the system eventually becomes 
an anisotropic = 1 IQH state. However, for intermedi- 
ate 9 value, we find that crystal-like order develops in the 
density profile of the system, suggesting the competition 
between various parameters and phases. By decompos- 
ing the dipolar interaction into isotropic and anisotropic 
components, we provide a simple explanation to quan- 
titatively understand the phase transitions induced by 
anisotropy. The various competing orders and ground 
states are summarized in a complete phase diagram in 
the parameter space spanned by the tilt angle and the 
confinement. 

In the presence of anisotropy in the interparticle inter- 
action, we lose the rotational symmetry when treating 
the system in a disk geometry, hence the total angular 
momentum is not a good quantum number any more. 
Nevertheless, we demonstrated that one can still compute 
the expectation value of the total angular momentum op- 
erator for eigenstates and use it to characterize the vari- 
ous ground states emerged as the results of competitions 
between confinement and anisotropy, as well as between 
the isotropic and anisotropic components of the interpar- 
ticle interaction. Particularly interesting is that the re- 
silient features of topological order, such as the presence 
of hierachical FQH ground states and the low-energy ex- 
citations pertaining to the density deformation along the 
edge of an incompressible quantum Hall droplet, remain 
robust in the presence of anisotropy. While we showed 
that the calculation of the fidelity of the ground state 
can be used as a probe to detect phase transitions be- 
tween states with different topological order, we believe 
the mean angular momentum treatment can be readily 
generalized to the calculation of, e.g., entanglement spec- 
trum , which also facilitates the detection of topolog- 
ical order. 

In this paper we demonstrated that an incompress- 
ible FQH state with a large excitation gap can survive 
a fairly large amount of anisotropy. This implies that 
the Laughlin state, the exact ground state produced by 
the hard-core potential, is stable against the anisotropic 
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perturbation. The finding is not unexpected given the 
remarkable stabihty of the Laughhn states in the pres- 
ence of long-range interaction, finite thickness of the two- 
dimensional electron gas, and disorder [29*1. It is, how- 
ever, intriguing to ask the effects of anisotropy on more 
exotic non-Abelian quantum Hall states. A more chal- 
lenging question would be whether it is possible, by tun- 
ing the anisotropy, to enhance a certain FQH state, hope- 
fully with exotic statistics. The present paper paves a 
path toward these questions. 
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